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ABSTRACT 

This is the second paper in a series reporting on results from a survey conducted with the 
ESO VLT/X-shooter spectrograph. We target high metallicity damped Lyman-a absorbers 
(DLAs) with the aim of investigating the relation between galaxies detected in emission and 
those detected in absorption. Here, we report on the discovery of the galaxy counterpart of the 
z abs = 2.58 DLA on the line-of-sight to the z = 3.07 quasar SDSS J 091826. 16+163609.0 
(hereafter Q 0918+1636). The galaxy counterpart of the DLA is detected in the [Olll] A5007 
and [Oil] AA3726, 3729 emission lines redshifted into the NIR at an impact parameter of 
2.0 arcsec (16 kpc at z = 2.58). Lya emission is not detected down to a 3u detection limit 
of 5 x 10~ 18 erg s _1 cm~ 2 , which, compared to the strength of the oxygen lines, implies 
that Lya emission from this galaxy is suppressed by more than an order of magnitude. The 
DLA has one of the highest metallicities measured so far at comparable redshifts. We find 
evidence for substantial depletion of refractory elements onto dust grains. Fitting the main 
metal line component of the DLA, which is located at z a b s = 2.5832, we measure metal 
abundances from Znll, Sll, Sill, Cm, Mnll, Fell and Nill of -0.12 ± 0.05, -0.26 ± 0.05, 
-0.46 ± 0.05, -0.88 ± 0.05, -0.92 ± 0.05, -1.03 ± 0.05 and -0.78 ± 0.05, respectively. In 
addition, we detect absorption in the Lyman and Werner bands of molecular hydrogen (H2), 
which represents the first detection of H2 molecules with X-shooter. The background quasar 
Q0918+1636 is amongst the reddest QSOs at redshifts 3.02 < z < 3.12 from the SDSS 
catalogue. Its UV to NIR spectrum is well fitted by a composite QSO spectrum reddened 
by SMC/LMC-like extinction curves at z a b s = 2.58 with a significant amount of extinction 
given by Ay « 0.2 mag. This supports previous claims that there may be more metal-rich 
DLAs missing from current samples due to dust reddening of the background QSOs. The fact 
that there is evidence for dust both in the central emitting regions of the galaxy (as evidenced 
by the lack of Lya emission) and at an impact parameter of 16 kpc (as probed by the DLA) 
suggests that dust is widespread in this system. 
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1 INTRODUCTION 

The comparison between absorption-line selected and emission- 
selected galaxies at redshifts z > 2 has a long history (e.g., 
ISmith et al.|[l98^ llvfoller & Warred [l993l IWolfe et alj|2005t and 
references therein). However, progress in this field has been slow 
and, for many years, ther e has been l ittle overlap between ob- 
servational samples (e.g ., Fynbo et al. 19991: llvfoller et all |2002| ; 
IColbert & Malkanl 120021 ; iKulkarni et al.1 120061) . Recently, some 
progress has been made though in order to build the bridge between 
the two populations. Emission-selected galaxies have been studied 
to much deeper rest f r ame flux limits than a decade before (e.g., 
Sawicki & Thompson! 2004 iGronwall et al.l |2007|; iRauch et all 



2QQ8l:l0uchi et al.ll2008l;lGrove et al.ll2009l : lReddv & Steidelll2009l ; 
Cassata et alj|201 lh . iFvnbo et al.l J2010I hereafter F10) also pre- 
sented the first results from a small survey to search for the 
gal axy counterparts o f metal-rich damped Lya absorbers (DLAs; 
see I Wolfe et alj|2005l . for a review). In F10, the detection of the 
galaxy counterpart of a high-metallicity DLA at z a b s = 2.354 to- 
ward Q 2222— 0946 was presented. F10 described the strategy and 
sample selection of the survey in details. Here, we reiterate that 
the candidates were selected amongst SDSS QSOs based on the 
strengths of Sill and Fell absorption lines (i.e., not directly from 
Lya). 

In this paper, we present new results based on observations 
of the second target of the survey, the z — 3.07 quasar SDSS 
1091826.16+163609.0 (hereafter Q0918+1636). Q0918+1636 
was selected as its spectrum features a metal-rich absorber at 
z — 2.412 having rest frame equivalent widths (EWs) of the 
targeted SillA1526 and FeIlAAA2344, 2374, 2382 lines of 2.1 A, 
2.1 A, 1.4 A and 2.6 A, respectively. The properties of this ab- 
sorber will be the subject of another paper describing a sample 
of objects. For now, we note that this absorber, which is also a 
DLA, has a high metallicity with [Si/H] = —0.6. Here, we de- 
scribe the serendipitous discovery and properties of a second DLA, 
at z abs = 2.5832, along the line-of-sight to Q0918+1636. This 
system happens to also fulfill our sample selection criteria, with 
rest frame EWs of 2.4 A, 2.4 A, 1.6 A and 3.0 A for SillA1526 
and FeIlAAA2344, 2374, 2382, respectively. The Fell lines from 
this system are redshifted outside of the wavelength interval suited 
for automatic searches of absorption lines in SDSS spectra and 
hence were not identified. Interestingly also, neither of the two 
DLAs in this QSO spectrum have been i dentified in automatic 
DLA (i.e., Lyq-based) s earches in SDSS (Prochask a et alj|2005l : 
Noterdaeme et al. 2009b). This is due to a recently identified bias 
against the detection of damped Lyman-alpha absorptio n lines at 
the blue end of the spectra (see lNoterdaeme et al . 2009b). Indeed, 
the signal-to-noise ratio is decreased significantly by the presence 
of any strong and/or several closeby DLA lines. 

Throughout this paper, we assume a flat cosmology with 



* Based on observations carried out under prog. ID 084.A-0303 with the 
X-shooter spectrograph installed at the Cassegrain focus of the Very Large 
Telescope (VLT), Unit 2 - Kueyen, operated by the European Southern Ob- 
servatory (ESO) on Cerro Paranal, Chile, 
f E-mail: jfynbo@dark-cosmology.dk 



f^A = 0.70, f2 m = 0.30 and a Hubble constant of Ho = 
70kms~ 1 Mpc~\ 



2 OBSERVATIONS AND DATA REDUCTION 

Q0918+1636 was observed on February 16 2010 with X-shooter 
at the VLT. The QSO was observed at three position angles (PAs), 
namely 60°, —60° and 0° (all East of North). The purpose of using 
three slit positions is to cover a larger field of view around the QSO 
as shown in fig. 1 of F10. The integration time at each PA was 3600 
s, and a 1.3"-wide slit in the UVB arm and 1.2"-wide slits in the 
VIS and NIR arms were used. 

The expected resolving power with the above setup is 4000, 
6700 and 4300 in the UVB, VIS and NIR arms, respectively. This is 
confirmed by the width of sky emission lines in the spectra. How- 
ever, the seeing during the observations was significantly smaller 
than the widths of the slits (i.e., 0.67 arcsec in all three exposures 
as measured from the width of the QSO trace around 7500 A) and 
hence the true spectral resolution is higher than that measured from 
sky emission lines. In the VIS arm, we measure the resolution di- 
rectly from the width of telluric absorption lines. We then assume 
that the ratio between expected and true resolutions is the same 
in all three spectroscopic arms, and the degradation of seeing as a 
function of decreasing wavelength is calculated as the ratio of the 
wavelengths to the power of 0.2. In this way, we infer resolving 
power values of 6400, 11900 and 8800 in the UVB, VIS and NIR 
arm, respectively. 

We pro cessed the spectra us ing the X-shooter data reduction 
pipeline (see lGoldoni et alj |2006. see also F10). First, raw frames 
were corrected for the bias level (UVB and VIS) or dark current 
(NIR). Then, after background subtraction, cosmic ray hits were 
detect ed and removed using the method developed bv lvan Do kkum 
< l200ll) while sky em i ssion lines were subtracted using the method 
described in iKelsonl J200I) . After division by a master flat-field, 
the spectral orders were extracted and rectified in wavelength space 
using a wavelength solution previously obtained from calibration 
frames. The orders were then merged and in their overlapping parts 
the merging was weighted by the corresponding errors which were 
propagated in the process. From the resulting merged 2D spectrum, 
a one-dimensional spectrum of the QSO was extracted. This ID 
spectrum together with its error file and bad-pixel map are the final 
products of the reduction. Intermediate products such as the sky 
spectrum and individual echelle orders (with errors and bad-pixel 
maps) were also produced. The spectra were flux-calibrated using 
a spectrophotometric standard star observed during the same night. 
The flux calibration was checked against the flux-calibrated SDSS 
QSO spectrum and found to be consistent with it. 



3 RESULTS 

3.1 Emission properties of the DLA-galaxy counterpart 

We do not detect Lya emission at any of the three position angles 
(see Fig. [TJ. We calculate a conservative 3<j detection limit from 
a 1000 kms -1 x 2.0 arcsec square aperture inside the trough of 



the z = 2.58 damped Lya absorption line. The velocity width of 
1000 kms" 1 is based on the Lya line from the galaxy counterpart 
of the DLA towards Q 2222 -0946 (F10, see also theoretical pro- 
files in lLaursen et al.l d2010h ). whereas the spatial width of 2 arcsec 
is based on the fact tha t Lya emission can be quite extended (e.g., 
lM0ller & Warrenll 19981 ; iFvnbo et al.ll2003l ; lRauch et al.ll2008h . The 
resulting limit is 5 x 10 -18 erg s -1 cm" 2 . Assuming that all Lya 
photons escape, case-B recombination and the rel ation between Ha 
luminosity and star formation rate (SFR) from iKennicuttl dl998h . 
we would place an upper limit on the SFR of 0.3 Mq yr" 1 . Un- 
fortunately, the expected position of the Ha emission line is in a 
part of the NIR spectrum where the sky background is too high 
to allow for a useful detection limit to be derived. However, we 
do detect the [Om] A5007 emission line in the PA= 60° spec- 
trum at an impact parameter of 2.0 arcsec (see Fig.QJ. The redshift 
of the [Om] A5007 line is consistent with that of low-ionisation 
metal lines to within an uncertainty of about 50 kms" 1 . This line 
is detected neither in the PA= —60° nor the PA= 0° spectra. We 
also detect both components of the [Oil] AA3726, 3729 doublet but 
with a lower signal-to-noise ratio and in a region somewhat affected 
by telluric absorption lines. The flux of the [Om] A5007 line is 
/ = 1.7 ± 0.2 x 10" 17 erg s -1 cm -2 , which corresponds to a 
luminosity of L — 9 x 10 41 erg s" 1 . This is a lower limit due to 
the possibility of slit-loss. For the [Oil] doublet, we infer a total 
flux of about 2.5 x 10" 17 erg s _1 c m~ 2 after correctin g for telluric 
absorption. Using the relation from IKennicuttl |l998), we infer a 
SFR of about 20 Mq yr" . Hence, Lya emission from this system 
appears to be suppressed by more than an order of magnitude. 

In order to exclude the possibility of the existence of a galaxy 
counterpart located at small impact parameter, we also performed 
spectral point spread function subtraction as in F10. No emission 
line is detected at an impact parameter smaller than 2 arcsec. Any 
line with a flux similar to or larger than the detected [Om] line 
would be detected in the data. 



3.2 Absorption-line properties of the DLA 

3.2.1 Metal lines 

The total Hi column density of the system is well constrained by the 
damped Lya absorption line to be log JV(Hl) = 20.96 ± 0.05 (see 
Fig. [TJ, where the error is a conservative error on the uncertainty, 
which is dominated by the systematic error from the normalisation. 
In particular, the fact that the DLA line is located in the red wing 
of the Ovi emission line of the QSO complicates the normalisation 
somewhat. 

The absorption profiles of the strongest Fell and Mgll lines 
span a velocity range of up to 600 kms -1 (see FeIlA2382 in the 
upper right panel of Fig. |2j- However, most of the absorption as 
seen from unsaturated or weakly saturated lines from singly ion- 
ized species is concentrated within a single broad velocity compo- 
nent at z a bs = 2.5832(2). In order to derive overall metallicities in 
the neutral gas phase, we focus on fitting this main component. For 
the fitting, we use the package FITLYMAN as available in MIDAS. 
From the weak SillA1808 line, one can see that this main compo- 
nent is responsible for 85% of the total EW and hence most likely - 
at least - 85% of the total column density of low-ionisation species. 

The results of Voigt-profile fitting are summarised in Table Q] 
and displayed in Fig. [2] The instrumental resolution FWHM is 47 
and 25 kms -1 for the UVB- and VIS-arm spectra, respectively. 
With a broadening parameter b — 52 kms -1 , the fitted metal-line 
profile is well-resolved. However, should it include in reality one or 
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Figure 1. Top panels: The detected [Oil] AA3726, 3729 and [Om] A5007 
emission lines from the DLA-galaxy counterpart observed 2.0 arcsec away 
from the QSO continuum. The arrows point to the predicted centroids of 
the lines if these were at the same redshift as the low-ionisation absorption 
lines from the DLA. Middle and bottom panels: The region around Lya 
from both the two- and one-dimensional spectra with PA= 60°. As can be 
seen, the Lya emission line is not detected. The red dashed and dotted lines 
show the Voigt profile fit to the damped Lyce line and associated 1<t uncer- 
tainties, respectively, corresponding to a neutral hydrogen column density 
of log JV (Hi) = 20.96 ± 0.05. 



several narrow components, hidden saturation could be an issue and 
strictly speaking our measurements should be considered as lower 
limits. This means that the overall metallicity of this DLA could be 
even higher than solar. However, apart from SlI and to some extent 
also Sill, we are considering relatively weak absorption lines in the 
fitting process (see Fig. |2]l so that the measured column densities 
for the corresponding species should depend weakly on individual 
b values if at all. 

The usual complications of line blending in low-resolution 
spectra are avoided here as a single well-defined component is fit 
and therefore the presence of blending with unrelated absorption 
line features is easier to identify. Absorption from neutral carbon 
is detected (QAA1560, 1656) but due to blending with Cl* and 
Cl** lines neutral carbon is not considered in the analysis any fur- 
ther. Also, it is not possible to infer whether ClI* is present due to 
blending with the CIIA1334 line which is highly saturated. Hence, 
we cannot determine the relative strengths of ClI* and ClI and use 
them as diagnosti cs to infer the presence of local radiation fields 
dWolfe et alj|20Q8h . The errors given in Table Q] are the formal sta- 
tistical errors from FITLYMAN. The systematic errors from, e.g., 
normalisation and not knowing the instrumental resolution pre- 
cisely, are larger than this. An error of 0.05 dex would be a con- 
servative estimate including all error sources. 

The metallicity of the system is exce ptionnally high. From 
Znll, which is little depleted onto dust grains dMever & RothTl 990), 
we infer a metallicity of —0.12 ± 0.05. Note that the splitting be- 
tween ZnIlA2026 and MglA2026 is 50 kms" 1 implying that Mgl 
cannot significantly contribute to the Znll column density. Here, 
we ad opt the solar photosphere abundances from lAsplund et al.l 
(2009). Given that we only fitted the main velocity component, 
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the overall metallicity of the DLA must be solar or higher than 
solar in case of hidden saturation. We also note that the [Mn/Fe] 
ratio is lar ger than solar, whic h is only seen at or above solar 
metallicity dLedoux et alj | 2002|). For elements which are sensi- 
tive to dust depletio n (e.g jMever & Rofhlll99fj| : |Pettini et aljl997l ; 
iLedoux et al.ll2002h . we find [Fe/H] = -1.03, [Ni/H] = -0.78, 
[Mn/H] = -0.92, and [Cr/H] = -0.88 implying substantial dust 
depletion for these elements. 

To compare the kinematics of the absorption-line profiles with 
that o f other DLAs, we follow the procedure of ILedoux et al.l 
(2006) and calculate the line-profile velocity width, AV, as 
c [A(95%) - A(5%)] /%, where A(5%) and A(95%) are the wave- 
lengths corresponding to, respectively, the 5 and 95 percentiles of 
the apparent optical depth distribution, and A p is the first momen t 
(the average) of this distribution (see fig. 1 of ILedoux et alj|2006l) . 
We again choose the Sill A 1808 transition as it is a low-ionisation 
transition and the line is only mildly saturated. The apparent line 
optical depth and the derived velocity width are shown in Fig. [3] 
We infer a velocity width of 295 km s~ in good agreement with 
the velocity-metallicity relation for DLAs ILedou x et alj|2006l) . 



3.2.2 Molecular absorption 

Absorption lines from the Lyman and Werner bands of molecular 
hydrogen (H2) are detected at z — 2.5832 in the X-shooter spec- 
trum of Q0918+1636 (see Fig.[4]l. The observed velocity extent of 
the H2 profile is wider than the instrumental resolution FWHM (48 
km s" 1 at 3800 A) and much wider than the typical Doppler param- 
eter of individual molecular lines as m easured from high-r esolution 
spectroscopy (usually b ~ 3 kms^i fLedoux et alj|2003l) . This in- 
dicates that the actual H2 profile comprises several blended com- 
ponents. Two absorption peaks separated by about 55 kms -1 can 
be distinguished in the H2 profile, while they are not seen in the 
profile of metal lines that are fitted with a single component having 
b = 52.1 kms -1 . This means that the metals are somehow continu- 
ously spread over FWHM = b x 2^/(ln 2) ~ 87 km s _1 , while the 
H2 components are discrete and well-separated in velocity space. 
Indeed, this is a general behaviour as seen in other H2-bearing sy s- 
tems studied at high spectral resolution (e.g. ILedoux et al . 2002b) 

In order to estimate the H2 column density, we compared the 
data with synthetic profiles of H2 absorption consisting of two com- 
ponents separated by 55 km s _1 built with the same spectral resolu- 
tion and binning. We varied the column density in rotational levels 
from J = to J = 3 for each component until we got a synthetic 
spectrum consistent with the data. We repeated this exercise for a 
wide range of Doppler parameters (b — 1 to 20 km s" 1 ) to get a 
realistic range on the total H2 column density: iV(H2)=1.5x 10 16 - 
1.1 x 10 19 cm -2 . This corresponds to an overall molecular fraction 
/ = (2iV(H 2 )/(2iV(H2)+iV(Hi))) = 3.3 x 10~ 5 -2.4 x 10~ 2 . 

Details on the column densities in each rotational level and 
for the two components are given in Table [2] These measurements 
are only indicative and should be considered with caution because 
of the various uncertainties. Only high spectral resolution obser- 
vations are suitable to properly resolve the H2 line profile in indi- 
vidual components, measure accurately the column densities, and 
study the physical cond itions in the gas (e.g. JSrianand et all 20051 ; 
Noterdaeme et al. 2007). Nevertheless, the presence of H2 in this 
system is firmly established. 
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Figure 2. Results of Voigt-profile fits to low-ionisation lines from the 
z a bs = 2.5832 DLA toward Q0918+1636. Most of the absorption is con- 
centrated within one velocity component. In order to derive overall metal- 
licities in the neutral gas phase, we focus on fitting this main component 
which is responsible for at least 85% of the total metal column density (see 
text). The Voigt-profile fits are shown with a red line and the residuals from 
the fits are shown with a blue line (for display purposes, a constant of 0.5 
has been added to the residuals in some of the panels). 
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Table 1. Ionic column densities in the main metal-line component of the 
2 aba = 2.5832 DLA toward Q0918+1636. 
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Figure 3. The profile of the Sill A1808 line. The profile consists of a main, 
broad component and a wing towards the red end of the profile. The velocity 
width of the line, measured following the method of Ledoux et al. (2006), 
is 295 kms - 1 (marked with a horizontal segment). Ai and A2 are the start 
and end wavelengths used to integrate the profile. 



Table 2. H2 column densities. The velocities, Av, of the two components 
are given in kms" 1 with respect to 2 a b s = 2.5832. 
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Figure 4. Portion of the X-shooter spectrum featuring H2 absorption lines 
superimposed on the Lya forest. The label over each absorption line indi- 
cates the band (L: Lyman, W: Werner), the vibrational level of the upper 
state, the branch (P, Q or R for A J = +1, 0, —1, respectively) and the 
rotational level J of the lower state. The synthetic profiles corresponding to 
b = 1 and b = 20 kms -1 are overplotted as a shaded area. The bottom 
panel represents the stacking of unsaturated, unblended lines (L0R0, L1R0, 
L1R1, L4R0 and L8R0), revealing better the two peaks in the H2 velocity 
profile. 



4 DISCUSSION AND CONCLUSIONS 

4.1 Nature of the DLA galaxy counterpart 

The z abs = 2.5832 DLA toward Q0918+1636 is amongst 
the highest me tallicity DLAs kno wn to date (compared, e.g., to 
the sample of iKaplan et alJ I201LT) . The fact that we detect the 
galaxy counterpart of this DLA is consistent with our suggestion 
that h igh metallicity DLAs should have bright g alaxy counter- 
parts (Ma iler et al.ll2004l: iFvnbo et alJ I2008L 1201(1) . Indeed, stel- 
lar mass, metallicity and star formation rate in galaxies are in- 
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timately related ( Mann ucci etal] l201Ch . Galaxy mass-metallicity 
relati ons have been found to exist up to at l east z = 3.5 
(e.g.. iTremonti et al j 120041; ISavaglio et al] 120051 : lErb et alj|2006l: 
Kewlev & Ellisorj 120081 ; iMaiolino et al] 1200% lLamareil le et all 
20091) . The mass-metallicity relations at high redshift are steeper 



and offset towards higher stellar masses for a given metallic- 
ity compared to the relation obser ved in the local Universe (see, 
e.g., rig. s in Maio linoetali r2008). Therefore, by selecting high- 
metallicity systems we should pick the most massive structures. In 
turn, because the star formation rate in galaxi es is at any redshift 
correlated with total stellar mass (see fig. 1 in button et al]|201(A 
and references therein), the most massive galaxies have expected 
star formation rates high enough (SFR> 10 Mq yr" 1 ) to be de- 
tected with our observing method. It is encouraging that we have 
already found a second case after F10 in agreement with this rea- 
soning. 

The large impact parameter of 2.0 arcsec, corresponding to 
16.0 kpc, is also consistent with the expectations from a simple 
model (see fig. 3 in lFvnbo et al . 2008). The detection of H2 absorp- 
tion is consistent with the previousl y found trend th at H2 is predom- 
inantly found in high-metallicitv jPetitiean et alj|2006l) and dusty 
dLedoux et alj|2003t iNoterdaem e et alj|2008l) systems. Unlike the 
system studied in F10, Lya emission from the z a b s = 2.5832 DLA 
galaxy toward Q 0918+1636 is strongly suppressed. A plausible 
reason for this is dust obscuration given the high metallicity and 
the evidence for d ust in the galaxy alo ng the line-of-sight to the 
background QSO dCharlot & Fallll 19911) . Whereas we have no di- 
rect information on the amount of dust in the centre of the galaxy, 
presumably traced by the centroid of the line emission, it is likely 
to be higher than in the outer parts of the galaxy traced by the DLA. 



4.2 Lya escape fraction 



Whereas the limit we infer for the Lya luminosity is a simil ar to the 
mean Lya lumi nosity for DLA galaxies in ferred bv lRahmani et al] 
d2010l) (see also lRauch& Haehneh1 l20irjh it is still very low com- 
pared to the strength of the [Olll] and [Oil] emission lines from the 
galaxy. 

The escape fraction / esc of Lya photons at high redshift is 
a highly coveted quantity as a plethora of interesting physical 
properties of galaxies can be revealed by studying this line. A 
number of authors have attempted to assess this matter, yielding 
quite different results; the standard procedure has been to com- 
pare the ratio of Lya-inferred SFRs to those of other proxies with 
wha t is expected theoretically as suming case-B recombination, 
e.g., iLe Delliou etal] J200l |200d a universal 2% from compar- 
ing observational data with modelled galaxies), lDavaletalj|2009l, 
30% from comp aring simulated galax ies with observed luminosity 
functions), and Gawi ser et al] d2006l . 80% from comparing Lya- 
inferred SFRs to spectr al energy distribution ( SEP) modelling of 
observed Lya emitters). lLaursen et aj ]2009bh, us ing the Lya ra- 
diative transfer code MoCaLaTA jLaursen et al]|2009ah , calcu- 
lated the full radiative transfer equations in simulated galaxies, ob- 
taining an anticorrelation between / csc and galactic size, from order 
unity to order a few per cent. 

Similarly, althoug h lHaves et al] feoich found that an average 
/csc = 5.3% adequate to explain their observations, their fig. 3 
shows a trend of / csc decreasing with E(B — V). Comparing the 
inferred SFR from [On] with that of Lya, an upper limit of ~ 1% 
is inferred for the galaxy counterpart of the DLA studied here. The 
extinction is found to be Av — 0.2 (see below), corresponding 
to a colour excess of E(B — V) ~ 0.06. The extinction in the 



centre of the gala xy is likely to be higher. Furthermore, in calculat- 
ing E(B — V) lHaves et al] ( |2010|) assume a metallicity of 1/3 so- 
lar, whereas the present galactic system was shown to have at least 
solar metallicity, thus probably resulting in higher extinction. An- 
other possible cause could be that SFRs inferred from [On] are after 
all associated with rather large uncertainties, as the mean [Oll]/Ha 
ratios vary substantially from ga l axy to galaxy, u p to an order of 
magnitude l lGallagheretaiIl9 89; Kenni cuttlll992l) . 



4.3 Dust reddening 

In Fig. [5] we show the X-shooter spectrum of Q0918+1636. 
With a dashed lin e we show the composite QSO spectrum from 
iTelfer et al] d2002l) and it is evident that Q0918+1636 is redder 
than the composite spectrum. The evidence for dust in the sys- 
tem comes from the observed strong depletion of refractory el- 
ements. In addition, the spectrum of Q0918+1636 appears sub- 
stantially redder than the spectra of typical QSOs around z = 3. 
Q0918+1636 was selected as a candidate high-z QSO based on 
its red it* — g* colour. It falls outs ide the selection crite ria for can- 
didate z < 3.6 QSOs in SDSS l lRichards et al]|2003h . In Fig. [6] 
we compare the colours of Q0918+1636 with the colours of QSOs 
with redshift within 3.07 ± 0.05 from the SDSS DR7 QSO cata- 
logue after removing BAL QSOs ( Schneider et al. 2007). As seen, 
the distributions of QSO colours are relatively narrow. However, 
the distributions have small extensions towards redder colours and 
Q0918+1636 is among the reddest 4%, 0.7% and 12% of QSOs 
at these redshifts for u* —g* , g*—r* and r*—i*, respectively. With 
contours we show t he colours of the stellar locus from stripe 82 
in the same colours Uvezic e t al. 2007). Only in the u*—g* colour 
has the reddening moved Q 0918+1636 further away from the stel- 
lar locus and in these bands the QSO is very faint. Hence, it is 
very difficult to select reddened QSOs at these redshifts using op- 
tical colours. In a future study, we will address how one can se- 
lect even redder QSOs using photometry in the near-IR bands. The 
most likely cause of the red colour of Q 0918+1636 is dust in the 
z a bs = 2.5832 DLA given its high metallicity and strong evidence 
for depletion of refractory elements on dust grains. Dust in the 
2 = 2.412 DLA could also contribute, but as mentioned in the in- 
troduction this DLA not as metal rich as the z a b s = 2.5832 DLA. 
The presence of dust along the line-of-sight is supported by the fact 
that the SED of Q 09 1 8+ 1 636 can be w ell fitted with the composite 
QSO spectrum from ITelfer etal] J2002l) reddened by SMC/LMC- 
like extinction at z a b s = 2.5832 with a total amount of extinction 
gi ven by Ay = 0.2 mag. Here, we have used the prescription given 
in |Pe3(l992) to model the extinction curve. The agreement with the 
observed spectrum is very good. It is not possible to get a good fit 
with the MW extinction curve. 

The corresponding extinction-to-gas ratio, Av /AT(Hl) ~ 
2.5 x 10 -22 mag cm 2 , is i n between the rati os inferred from 
LMC and MW sightlines (e.g. lGordon et alj2003l) and substantially 
higher than what i s found on average for DLAs (A v /N(Hl) ~ 2- 



4 x 10~ 23 mag cm 2 rVladilo et alj2008h . The fact that reddening has 
moved Q0918+1636 to the boundary of the colour space where 
QSOs are selected in SDSS is consistent with previous studies ar- 
guing that metal-rich DLAs most likely are systematically miss- 
ing from curre nt DLA samples from o ptical QSO surveys (e.g. , 



Pei et al. 1999 and references therein; Noterd aeme et al] 2009a 
20ld ; |Pontzen & Pettinill2009h . 



In conclusion, our results support the suggestion of F10 that 
high-metallicity DLAs are associated with bright galaxy counter- 
parts. The 2 a bs = 2.5832 DLA has a high metallicity (for its red- 
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Figure 6. Here, we compare the colours of Q0918+1636 with the colours 
of other SDSS QSOs within Az = 0.05 of the redshift of Q0918+1636 
(BAL QSOs excluded). The top two panels show colour-colour diagrams 
and the lower three panels show the histograms of colours. The con- 
tours in the upper panels show the colour distribution of the stellar locus. 
Q0918+1636 (marked as a circle in the colour-colour plots and a dashed 
line in the histograms) is among the reddest of all SDSS QSOs in this red- 
shift range (4% in u*—g*, 0.7% in g*—r* and 12% in r*—i*). 



shift) of close to solar despite probing a region about 16 kpc away 
from the central emitting region of the DLA galaxy counterpart. We 
establish the presence of H2 absorption, and also dust grains based 
on the depletion of refractory elements and from reddening of the 
background QSO. 
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Figur e 5. The sp ectrum of Q 0918+1636 after flux calibration. The overall shape of the spectrum is well fitted by the composite QSO spectrum from 
iTelfer et all (2002). In the figure, the unreddened composite spectrum is shown with a dashed line, and the same spectrum reddened by SMC- and LMC-like 
extinction curves with rest frame Ay = 0.2 mag is shown with solid red and blue lines, respectively. The inferred dust-t o-gas ratio is between those of the 
LMC and the MW. The spectrum has been corrected for galactic extinction with E(B — V ) = 0.025 from lSchlegel et alj ll998l) . 
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